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Abstract 

A summary of the successes of and obstacles to the gauge technique (a non- 
perturbative method of solving Dyson- Schwinger equations in gauge theories) is 
given, as well as an outline of how progress may be achieved in this field. 



1 The gauge technique 

Nature seems inordinately fond of gauge theories. Chromodynamics, elec- 
troweak theory and gravitation are all based on the gauge principle, featuring 
the groups SU(3), SU(2)xU(l) and SL(2,c), respectively. It is quite likely 
that any unified gauge model of the fundamental forces will also be a gauge 
field theory and it is conceivable that a supersymmetric version, in some 
higher dimension, will be founded on a local supergauge principle, although 
nature seems reticent to display her supersymmetric hand in the low energy 
spectrum of states. 

Whenever a field theory is invariant under local group transformations, 
the resulting Green functions obey a series of 'gauge covariance' relations 
which specify how the functions change under a variation of gauge]]. In QED 
these relations were originally found by Landau and Khalatnikov (1956, LK 

1 Purely photonic processes and amplitudes involving on-shell fermions remain gauge- 
invariant of course. 



1 



for short); later on, Johnson and Zumino (1959) and Zumino (1960) rederived 
them using functional methods. An immediate consequence of these ideas 
is the famous divergence properties of Green functions at each vector leg, 
named Ward- Green- Takahashi (WGT) identities after their early discoverers 
(Ward 1950; Green 1953; Takahashi, 1957). 

Heroic efforts have been expended in determining the behaviour of off- 
shell Green functions. This program has be carried out to high accuracy 
in perturbation theory for the standard model of particle physics and the 
results (at least for asymptotic momenta in QCD) are quite reliable because 
the perturbation series 'converge' to the order in which they have been cal- 
culated so far and also because the model is renormalizablef]. An alternative 
approach is to solve the full set of Dyson-Schwinger (DS) equations connect- 
ing the various Green functions by some suitable truncation procedure. For 
instance, one simple way is to sum the 'rainbow diagrams' of the self-energy 
or the 'ladder diagrams' of the vertex, which amounts to simplifying the DS 
equations for the two- and three-point functions in a particular way. By this 
means one can already discern the phase structure of a model as the coupling 
runs through certain critical values. 

In gauge theories, thanks to the existence of WGT identities between 
successive n-point Green functions, a more sophisticated approach suggests 
itself. The identities allow one to determine (longitudinal) pieces of the Green 
function in terms of the function with one less vector particle^. Hence if one 
discards transverse parts of the Green function (orthogonal to the photon 
momentum) or at least ties them in a particular fashion to the longitudinal 
parts, then this amounts to a truncation which produces a closed set of 
equations among the functions — to some level of approximation depending 
on n. For want of better terminology, we will coin this procedure the 'gauge 
technique' after Salam (Salam 1963; Salam and Delbourgo 1964; Strathdee 
1964) — even though it was originally tied to a specific way of 'solving' the 
gauge identities and began at level n — 2. Anyway, armed with solutions of 

2 They will eventually go wrong because the series in e 2 is believed to be only an 
asymptotic one. 

3 There is even a 'curly' version of such identities, obtained by Takahashi; they are 
rather more complicated than the ordinary divergence and we will have little to say about 
them in this paper. However in 2-dimensions, it should be noted that they are equivalent 
to the axial gauge identities (Delbourgo and Thompson 1982; Thompson 1983; Kondo 
1997). 
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the DS equations, one may go on to calculate quantities of physical interest, 
such as decay constants of hadrons (Maris and Roberts 1997; Ivanov et al 
1998). This approach (Thompson and Zhang 1987; Roberts and Williams 
1994) can be carried out in any number of space-time dimensions D and for 
any physical fermion mass m and it is quite interesting to see how the results 
depend on D and m separately, in various background (Cornwall 1986; Maris 
and Roberts 1998) configurations (temperature, field, etc.). 

In this paper we shall focus primarily on the archetype gauge theory, 
QED. This has the advantage of yielding simple abelian identities and of 
bringing into relief the principal obstacles which the gauge technique must 
circumvent before it can be considered a fully-fledged method, with results 
that are above reproach. The purpose of this paper is to highlight some of 
these problems and to suggest possible ways in which some progress may be 
achieved. 

Since we are dealing with gauge theories, the first issue we face is the 
gauge covariance of the technique, i.e. how well do the Green functions solu- 
tions obey the LK relations. The second question concerns the comparison 
with perturbation theory (in the fine structure constant a), i.e. how closely 
do the non-perturbative answers coincide with the perturbation series if we 
try to expand them in a. The third matter refers to the level of approxima- 
tion, i.e. to what extent do the results change with the level of truncation n. 
Since the technique concentrates primarily on the charged particle functions, 
the fourth issue is what happens to vector propagation, i.e. how will vac- 
uum polarization pan out in such schemes: this is an serious matter because 
photon renormalization, not electron renormalization, governs the momen- 
tum dependence of the running charge. Indeed vertex and charged particle 
renormalizations are largely irrelevant: they are gauge-dependent and can- 
cel out, leaving little physical trace. It means that going from quenched to 
unquenched solutions represents a major step physically. 

If these issues can be clarified in QED, there are definite lessons (Pa- 
pavassiliou and Cornwall 1986; Alekseev 1998; Cahill and Gunner 1998) for 
QCD and gravity. 
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2 Gauge identities and SD equations 



First we set up the framework for the ensuing discussion by summarising 
quickly the identities and the coupled equations for Green functions in QED. 
If nothing else this will fix our notation. Let G Xl}X2r (pi,p2, .; ki, &2, •) stand 
for a connected time-ordered Green function, with electron momenta p and 
photon momenta k, as labelled. It is usual to write G(p) = S(p) and G fMU (k) = 
D^ v {k) for the two-point electron and photon propagators respectively. It is 
convenient often to define the one-particle irreducible Green functions T(p\ k) 
by dropping all the pole parts of G(p; k) and by multiplying out by the 
inverse propagators of the external lines; for instance, the '3-vertex function' 
T{p',p; k) with k = p' —p arises via G x (p',p; k) = S(p f )T fl (p',p; k)S(p)D tlX (k), 
while the 'Compton part' T(p',p; k',k) of electron photon scattering occurs 
in 

G XK (p',p;k',k) = D Xu (k')S(p')[T^(p',p;k' } k)- 

T u (p',p' + jfc'; k')S{p + k)T^(p + k,p; k) - 
T v (p-k',p; k')S(p'-k)r^p',p'-k;k)}Sip)D„ K (k)(l) 

etc. The x-spa.ce Fourier transforms are obtained as convolutions; thus 

G\(x, y; z) = J d A x'dSj'd A z' 'D Xfl (z - z')S(x - x')F»(x' \y'; z')S{y'-y), 

and so on. The reason why we have mentioned the coordinate space version 
is because the gauge covariance identities, to which we shall presently turn, 
are best written in a;-space. 

The Dyson-Schwinger equations connect successive (renormalized) charged 
Green functions through the series^, 

S(p)(jf-mo) = Z~ 1 + ie 2 J d 4 kG x {p,p-k;k)-f X (2) 

G^(p',p;p' -p)(jf-mo) = S(p')>y^-ie 2 J d A kG ^ x (p , p - k; p' - p, k)^ x (3) 
etc. plus the uncharged cases starting with, 

Vllu = Z A D^(k) 0ll X D Xu (k) - te 2 Z Tt J d*pG^{p+k,p) lv . (4) 
4 Below and later on, we adopt the notation, d A k = d 4 k / (2ir) 4 , 5 4 (k) — (27r) 4 5 4 (fc), etc. 
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The nonlinearity of electrodynamics becomes evident via the infinite skeleton 
expansion of the higher-point Green functions, when one expresses the G in 
terms of their 1PI counterpart functions V. That is why one is obliged to 
invoke some kind of truncation for solving the equations nonperturbatively, 
instead of resorting to e 2 expansion — the idea being that a nonperturbative 
solution may reveal some dependence in 1/e 2 (usually in an exponent) which 
is not immediately apparent from the asymptotic e 2 series. 

All of the above applies, however one fixes the photon gauge. Now under 
a change of photon gauge, D^ v (x) — > D^(x) = D /Jil/ (x) — d^d u M(x), the 
renormalized Green functions are modified in a well-defined manner from G 
to G M . Thus 

S M (x) = exp[ie 2 M(x)}S(x), 

G^(x,y 1 z)=ex P [ i e 2 M(x-y)][G fl (x,y;z)+ i S(x-y)d;{M(x-z)-M(y-z)}], 

and so on. Let us call these the 'gauge covariance' or LK relations. It is 
simple to check that they are consistent with the SD equations in any gauge. 
A secondary consequence (though historically a primary one) is that the 
Green functions will satisfy the WGT identities, 

(p'-p) x G x (p',p;p - p) = S(p)- S(p), 

k^Gf, u {p',p] k', k) = G u (p',p + k; k') - G v {p - k,p; k'), 

whose soft /c-limit produces the Ward versions, G^(p,p;0) = —dS(p)/dp^, 
etc. The WGT identities also appear straightforward for the 1PI functions 
r, e.g. 

(p' -p) x T x (p',p;p -p) = S 1 (p)-S 1 (p). 

However the WGT identities are weaker than the LK relations, which can 
themselves become quite complicated for the amputated T, in contrast to 
those for the full Green functions G, especially when written in momentum 
space. It is therefore a vexing business to verify that any T, derived some- 
how in some M-gauge, properly obeys the required covariance identity; by 
comparison it is easier to investigate the matter for the non-amputated G. 
In fact if one only amputates the photon legs, the gauge covariance relations 
simplify a little further and make our task easier. For illustration, take the 
vertex function STS; one finds, like S M (x — y) = exp[ie 2 M(x — y)]S(x — y), 
that 

(STS) M (x,y,;z) =eM^ 2 M(x -y)}(STS)(x,y;z), (5) 
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and d z ■ (STS) M (x,y;z) = iS M {x - y)[5\x - z) - 8\y - z)\. (6) 
That completes the quick tour of the basis for the gauge technique. 

3 What gauge covariance implies 

Before tackling the spinor case, let us start with scalar electrodynamics, 
where the algebra and arguments are simpler. If A connotes the charged 
scalar propagator and V the fully amputated 3-point vertex in some gauge 
(specified by M = say), focus on the three relations: 

A M (x,y) = exp[ie 2 M(x - y)]A(x,y), 
(ATA) M (x,y;z) = exp[te 2 M(x - y)](AFA)(x, y;z), (7) 
d^AT^A) M (x,y;z) = iA M {x-y)[5\x-z)-5\y-z)l 

associated with the lowest functions in a different gauge M. Now in general, 
the off-shell 3-point function can be expressed in terms of two invariants, one 
associated with the longitudinal vertex and the other with a purely transverse 
vertex (k = p' — p): 

T,{p\p-k) = (p' + p)^( J 9'^/,A; 2 ) + [A;,(p'V)-(p / +p)^ 2 ]^(p'^p 2 ,A; 2 ), (8) 

where L and T are symmetric scalar function under p 2 <-> p' 2 . It follows 
that, when p 2 = p' 2 , the transverse part can be effectively combined with 
the longitudinal piece; this case applies in particular when one goes on the 
meson mass shell. 

It is the T part which largely governs^ the meson 'form factor' because 
the scalar WGT identity tells us unambiguously that the longitudinal part L 
cannot depend on the invariant (p' — p) 2 in any gauge; for it is always true 
that 

L M (p' 2 , P 2 ) = (A- 1A V) " A-' M ( P ))/(p' 2 -p 2 ). (9) 

Like L, the transverse contribution T must also change with gauge function 
M, but in a subtler way than L and one where the (p — p') 2 dependence 

5 It is worth noting that T has no singularities in the triangular variable fc 4 +p 4 +p' 4 — 
2p 2 p' 2 — 2p 2 k 2 — 2p' 2 k 2 , but that if we set p 2 ^ p' 2 there do arise logarithmic divergences 
(Ball and Chiu 1980) in perturbation theory as k 2 — > 0. 
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cannot be so easily forgotten. For suppose that in some gauge, we were to 
ignore T and wrote at the very least: 



(Ar M A)(p',p; k) = A(p')(p + P%A(p)L. 
Then in another gauge, according to (7), the result would get modified to 

(Ar M A)(p',p; k) = J d A k'A(p'-k')(p + p'-2k'),A(p-k')LE M (k'), 

where 

E M (k) = J A exp[ie 2 M(x) + life • x\. (10) 

A transverse amplitude T in the off-shell vertex would be ineluctably en- 
trained via the momentum numerator of the integrand. For instance, start- 
ing with first order perturbation theory, a transverse Lorentz-covariant lurks 
within the expression 

(AT A) M (v' v k) - [dV (P + P'- 2k % EM ( k ') (u) 

Only in the limit e 2 M = 0, when E M (k') = S\k'), will such a transverse 
term disappear. This argument teaches us two things about ensuring gauge 
covariance for general M: (i) that it is perilous to neglect the transverse 
parts of Green functions in non-perturbative treatments, and (ii) that one 
cannot always disregard the dependence of the amplitude on the momentum 
of the vector leg, i.e. one cannot purely^ use functions oip 2 ,p ;2 . 

The discussion increases in substance, richness and delicacy for fermions. 
Instead of one transverse and one longitudinal part, the vertex contains four 
independent longitudinal pieces and eight transverse pieces off-shell. In a 
real tour de force, these pieces have been computed by Ball and Chiu (1980) 
to first order perturbation theory in the Fermi-Feynman gauge, and by Kizil- 
ersu, Pennington and Reenders (1995) in any gauge. Being off-shell, the 
answers are extremely involved and we shall content ourselves with making 

6 This applies with force to various 'improvements' or corrections to the longitudinal 
vertex, consistent with multiplicative renormalizability, that have been suggested (Curtis 
and Pennington, 1993; Haeri 1993). In this connection it is worth remembering that the 
covariant gauge a — 3 produces a zero first order in a correction to the scalar propagator 
or to all orders, in the infrared limit. 
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three remarks: (i) there exist transverse covariants now which survive the on- 
shell spinor limit, such as ia^ v k v ', that have important physical consequences; 
(ii) no longer can one combine transverse vertices with longitudinal ones for 
physical fermions[|; (iii) there is a lot more freedom in 'solving' the WGT 
identity for the vertex, with distinct methods, all deemed free of unphysical 
singularities and all independent of the square of the photon momentum, 
yielding vertices differing by specific transverse terms T. (See the Appendix 
for details.) The various T arise automatically in the spinor version of (8) 
for arbitrary M, so it is futile to discard them in a general gaugeQ unless one 
gives up on the idea of satisfying the LK relations — a big disappointment for 
a gauge theory 

4 Comparison with perturbation theory 

If one does succeed in obtaining 'acceptable' solutions (presumably with an 
implicit dependence on e 2 ) of the DS equations by the gauge technique, then 
there is the obvious question of how well the answers stack up against pertur- 
bation theory, when expanded in powers of e 2 . At the very least one would 
hope that they would correct up to first order in a, as they are indubitably 
correct to zeroth order — but that hardly constitutes progress! It is not much 
good having them agree with perturbation theory in one gauge (for some L 
and T) but being wrong in another gauge. However that is what will likely 
happen from the transformation property of the G unless the various parts 
carry precisely the correct dependence on M and, as stressed previously, one 
countenances some k 2 dependence in STS. And it is no good avouching that 
the off-shell dependence of the propagator on the momentum looks 'reason- 
able', with 'suitable behaviour' in the infrared or ultraviolet limits, since one 
can change the behaviour at will, just by choosing the gauge function M 
however one likes. 

To clarify these points, let us consider covariant gauges, parametrized by 

7 In fact because physical fermions satisfy free equations of motion, the transverse co- 
variants are no longer independent but can be transmogrified into one another. 

Nevertheless one should recognize the privileged position of the covariant (Landau) 
gauge a — 0, because the first order in a correction to the wave function renormalization 
vanishes identically, as do all rainbow modifications of the spinor propagator. See Section 
6. 
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a real number a: 



D; u (k) = -V^/k 2 + (l-a)k,k v /k 4 . 

(The values a = 0, 1, 3 define the Landau, Fermi-Feynman and Fried- Yennie 
gauges respectively.) Upon introducing a regularization scale, like the elec- 
tron mass m, the Fourier transform of M(k) = —a/k 4 is obtained as iM(x) = 
— a ln(— m 2 x 2 ) / 16tt 2 . Therefore the coordinate space Green functions G of 
two charged fields separated by x will be multiplied by the gauge factor 

exp[ie 2 M{x)} = {-m 2 x 2 )~ ae ; e = e 2 /167i 2 = a/4n, 

with transform, 



E M 



(k) = \ d 4 x exp(tp-x)(-m 2 x 2 )' ai = _ 167r ^(2 ae) ^^^y^ 
J /c 4 r(ae) 



Expanding in powers of e 2 , one can check that 

E M (k) = 5\k) - te 2 a/k 4 + 0(e 4 ), 

with the second term on the right corresponding to the change of covariant 
gauge parameter a to first order perturbation theory. 

The transverse term has a complicated analytic form which only simplifies 
in various asymptotic limits; the behaviour in the ultraviolet regime is one 
such limit and King (1983) and Haeri (1988) put it to good use in correcting 
multiplicative renormalizability that is jeopardised in simple gauge technique 
ansatze; however, their procedure fails in the soft photon limit and does not 
do justice to the analytic behaviour in momentum transfer. To do better, it 
is useful to look at the form of T in scalar electrodynamics say, to first order 
in a or M as above. An examination of the Feynman graph integral, shows 
that it can be expressed in the parametric form 

-.2 



( p > 2 - p 2 )T = ^ J dadpd-y 5(1-^-/3-7) 



m 2 (l — a) — (p' 2 j3 + p 2/ ~f)a — (p — p') 2 ^ 
where the parametric numerator (for a = 1 say) is given by 

N = 2(p' 2 -p 2 )a + (P- 7 )[p 2 (l - P) +p f2 (l - 7) + m 2 (3 - a)] 



(12) 
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303 - 7 )/37[(p' 2 - P> + (P - PW ~ 7)]- 



The important thing is that the integral produces a result which vanishes for 
p 2 = p' 2 . 

It is interesting to see whether the answer can be rewritten as some kind of 
dispersion relation. To that end, consider the self-energy first. The standard 
Feynman parametric form leads to the scalar integral 



By changing variable to W 2 = m 2 /a, the expression above can be easily 
converted into the more familiar form 



Turning to the proper vertex part, if we make the change of variables in the 
Feynman parametric integral, (3 + j = a, f3 — j = au and a = m 2 /W 2 , we 
obtain an expression like 



p ,2 (l + u)/2 + p 2 (l - u)/2 + (p' - p) 2 (l - u 2 ){\ - m 2 /W 2 )(l - u 2 )/A - W 2 , 



and we see that only near the fermion mass shell (W 2 = m 2 ) and in the soft 
photon limit can one disregard the dependence on the momentum transfer. 
(In the ultraviolet limit where the spinor momentum is a long way from the 
p 2 = m 2 one can make other sorts of approximations.) 

What this exercise demonstrates is that under no stretch of the imagina- 
tion can one invoke a transverse vertex which is purely a function of p 2 and 
p' 2 . Even an invokation like 



misses the point altogether, since it fails to include the correct analytic struc- 
ture for a magnetic form factor. Remember also that the covariance LK rela- 
tion under change of gauge will inevitably create such structure, if none were 
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for the transverse part. The denominator D takes the form 
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initially present. Therefore we insist that any substantive improvements to 
the gauge technique ought to include these kinds of effects and agree with 
first order perturbation theory at the very least. Ansatze of the type, 

Tl(p',p)= dW dup T (W)K^/D(p',p,u,W) 

Jin J — 1 

stand some chance of capturing the main feature of transverse corrections. 

5 Level of truncation 

In an effort to improve upon the three-point Green function it is good to 
proceed to the next level of the DS equations. This will relate the three and 
four-point Green functions, via a Bethe-Salpeter like equation — see Parker 
(1984), Delbourgo and Zhang (1984). Note that in order to make use of the 
higher-point WGT identities, it is sensible to write this equation with one 
of the spinor legs preferred rather than the photon leg as is normally done; 
for with that choice one can use the four-point WGT identity to relate to 
relate that amplitude to the three-point Green function and thereby arrive at 
a self-consistent equation for the three-point function. (This is in complete 
analogy to the favourite way of handling the propagator and vertex function 
together and we can think of it as the n = 3 level improvement of the gauge 
technique, in contrast to the conventional n = 2 level). 

What is more to the point is that, at this level, the longitudinal and trans- 
verse vertices are treated on an equal footing while the two-point function 
(propagator) is obtained secondarily through the divergence of the three- 
point function. In practice, though, the equation for the vertex is very dif- 
ficult to solve, for one is not entirely sure of its full analytic representation, 
except through the perturbation expansion. The best one can do in such cir- 
cumstances is to use a double dispersion relation (in a Feynman parametric 
form say) and try to determine the spectral function self-consistently, but 
even that is not straightforward. The only progress to date has been the de- 
termination (Delbourgo and Zhang 1984) of the propagator spectral function 
in a manner which coincides with perturbation theory to order a 2 . Hopefully 
this problem will receive due attention in future. 
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6 Unquenching — vacuum polarization 



In most of the self-consistent calculations of the charged field propagator, 
the photon is taken as bare: this represents the so-called "quenched ap- 
proximation". Including the effect of vacuum polarization leads to a highly 
nonlinear equation for the propagator even in the best method of solving 
for the propagator, since vacuum polarization is determined by charged loop 
effects; that is one of the main reasons why the problem is normally avoided. 
The best that one can do in that situation is to make an inspired guess at 
the behaviour of the dressed vector propagator D(k) and (i) either include 
this in the computation of S(p), or (ii) introduce a running coupling in the 
interaction between the vector and the charged field. Now in QED the pho- 
ton receives small logarithmic corrections in the ultraviolet regime, but its 
low energy properties are largely unaffected; for that reason the quenched 
approximation is not too drastic a procedure for electrodynamics, at least in 
four dimensions^. But in QCD, the gluon exhibits asymptotic freedom in the 
ultraviolet — again a logarithmic correction — which hints at infrared slavery 
(confinement of colour?) at low energies and a propagator which is possibly 
more singular that the undressed form (Alekseev 1998; Cahill and Gunner 
1998). Thus a variety of models have been proposed and the corresponding 
quark S(p) found; the results depend critically on the assumed form of D{k) 
which is itself influenced by the gluon self-interactions and the ghost field 
effects. There is still some dispute about what is the correct form of D(k) 
and whether the quark propagator is an entire function of the momentum. 
We do not wish to get involved in these arguments; suffice it to say that the 
method with the best chance of being correct is the one that handles the 
gauge covariance properties correctly, in tandem with the gauge-invariant 
effective coupling g = g$Z{ Z 2 Z^' . 

In this connection it is worth recalling (Larin and Vermaseren 1993) the 
perturbative expansions of the beta functions for QED and QCD, which are 
of course gauge independent. Let e = g 2 /IQ-k 2 — ► a/47r for QED. Then 

oo 

/5(e) = e[ 73 (e) + 2 72 (e) - 2 7l (e)] = £ (3 n e n+l 

71=1 

In three dimensions charged loops alter substantially the low energy behaviour from 
D(k) = l/k 2 to D(k) oc ly/—k 2 , while in two dimensions the vector becomes massive 
through the Schwinger mechanism. 
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(In QED the anomalous dimensions are equal, 72 = 71, so (3 = €73 is deter- 
mined purely by the anomalous dimension of the photon field.) Up to three 
loops, 




325iV 2 
54 



5033iV 
18 




(13) 



where N stands for the number of charged fermions or flavours. The de- 
pendence on the gauge parameter a in QED arises through the anomalous 



Notice that the gauge dependence arises only at one-loop level and that it 
vanishes in the Landau gauge a = 0. More significantly there are higher 
order in a contributions which cannot be ignored; it is therefore fatuous to 
suppose that one can simply set the coefficient of fa in the spinor propagator 
equal to one in the Landau gauge — this is simply incorrect in higher orders. 

Anyhow it is reasonable to enquire what repercussions, if any, does the 
gauge technique have on vacuum polarization. There are two important 
points to check in this connection: whether the polarization tensor I\. liv {k) is 
transverse and whether it is gauge- independent. Since 



it is fairly clear that a non-perturbative calculation of II which dresses the 
fermions (S) but leaves the full vertex as bare (r — > 7), will not produce a 
transverse tensor; therefore this strategy is unacceptable. It is also easy to 
concoct a longitudinal approximation to the three-point Green function that 
does lead to transverse polarization, e.g. a mass-weighted spectral represen- 
tation, as described in the appendix. Carrying out the computations in the 
Landau gauge say, where 



10 See Larin and Vermaseren (1993) for a complete determination of the anomalous 
scaling functions in QCD, which are too long to reproduce for the purposes of this paper. 



scaling function for the spinorQ, 




LV(A0 = ie 2 Tr J d A p(ST^S)(p + k, V )% 
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(Sr L S)(x,y;z) = J dW p (W)S F (x - z\W)>yS F (z - y\W), (14) 



with S F (x I W) = (ry • d + PU)Air(a;|lU), we will obey gauge covariance by 
stipulating that, in any other gauge M, the expressions above are to be 
multiplied by exp[ie 2 M(x — y)\. This then leads to the gratifying result that 
the resulting vacuum polarization is transverse and gauge- independent, since 



if we approach the equal location limit x = y along a certain direction. This 
would imply that vacuum polarization is gauge independent and is given 
by the Landau gauge result, Hp V (z) = J dW p (W)Il llu (z\W), corresponding 
to a weighted mass integral. Actually this result is still not correct: the 
Green function (14) is insufficient to account for all higher order quantum 
corrections, because we must supplement it by a transverse (Landau gauge) 
contribution, which certainly does not spoil the transversality property. Un- 
less this is added, Pqed, will almost certainly be wrong. 

Thus, from all of this discussion, we see that the most pressing prob- 
lem in patching up the gauge technique is to incorporate transverse Green 
functions which have correct analytic and gauge- covariance properties. Until 
this is done, any physical results that are claimed to be a consequence of the 
technique are not to be fully trusted. 
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i,e 2 M{x-y) 



= Tr J dWp (W)S F (x - z\W)^S F (z - y\W) lv , 



(15) 
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Appendix 



In scalar electrodynamics, if A(p 2 ) = J dW 2 p(W 2 ) / (p 2 — W 2 ) stands for the 
meson propagator, the longitudinal expressions for the vertex, 

i>',p) = (P+^ A " 1( ^:^" 1(p2) , and 
(ArjA)(p', P ) 



are entirely equivalent. Although one may add any amount of transverse 
amplitude without affecting the WGT identity, it would be perverse to in- 
troduce such an additional piece unless there are good reasons to do so like 
reaching agreement with perturbation theory, eliminating sub divergences or 
trying to patch up the gauge covariance relation. This is precisely what 
has motivated King (in QED) and Haeri (in QCD) to incorporate particular 
transverse terms in the ultraviolet regime (King 1983; Haeri 1988). No such 
corrections are needed in the infrared regime, since all Green functions are 
effectively governed by the nonperturbative behaviour of the charged particle 
propagator (Delbourgo and West 1977a, 1977b; Delbourgo 1979; Atkinson 
and Slim 1979). 

In spinor electrodynamics the situation is much less clearcut because of 
the plasticity of the gamma-matrix algebra. Begin with the spinor propaga- 
tor, written in the equivalent forms, 

S(p) = J dWp(W)/(i>-W) or S'\p) =M(P 2 ) + B{p 2 ). 

For short, write A = A(p 2 ),B = B(p 2 ),A' = A(p' 2 ),B' = B(p' 2 ) and F = 
p 2 A 2 — B 2 , etc. Then there are at least three 'obvious' ways of 'solving the 
gauge identities', all of which are singularity free. From the proper vertex 
identity, one may factor out the momentum transfer and arrive at the first 
version (Ball and Chiu 1980), 

(B'-B)+ 1 -(fi>+ MA' -A) . 



A(p 2 ) - A(p' 2 ) 



= J dW 2 p(W 2 ) 



P 



/2 



w 2 



(p+p'] 



M p 2 - w 2 
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A second way is to carry out the factorization for the full Green function: 

~fl+j>{A' A\ ( B' B\ 



F> F \F> F 



A third way is to take advantage of the dispersive representation as a weighted 
mass integral and thereby arrive at the form 

These three versions are not identical to one another, in contrast to the scalar 
case. They differ from one another by particular transverse components 
(which of course have no effect on the WGT identity) and no version is more 
natural than any other at this level, unless other considerations intervene; 
thus they all behave smoothly as p 2 — > p 2 and they agree with lowest order 
perturbation theory. For instance, the difference between the first and second 
versions of the proper vertices can be expressed as 

A> -MM + B) + A% + U/A' + B')T, ( , A A + T,A 



(#A'-B') ^ ' "(M-B) 

where 2T M = 7^ — (fl— ^ l ){p'+p) fJ ,/{p' 2 —p 2 ) is a transverse Lorentz-covariant. 
Similarly, the third version can be rewritten in a more revealing form: 

(SriS)(j/,p) = J dWp(W) J ± w 2T„ J ± w + [S(p) _ s(p')] . 

Haeri and Haeri (1992) have shown that this particular spectral form of the 
longitudinal vertex can be converted into the equivalent but more elegant 
form 

which corresponds to the proper vertex solution 

j>(A' -A) + + 7m MB' ~B)+ 1 ^p' 2 A' - p 2 A) 



rj(p',p) = 



p' 2 — p 2 



featuring the inverse propagator functions A and B. 
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